clear; clc; close all;

u   = 0.1;
k   = 1.8;
a   = -1;
x0  = 0.1;
q0  = 0.1;
N   = 5000;    % 总迭代步数
discard = 100; % 舍弃前 500 步
keep    = N - discard;  % 用于绘制的步数

% 不同 b 值
bVals = [1.8, 18, 35];

for i = 1:3
    b = bVals(i);

    % 迭代系统, 返回 xArr, qArr
    [xArr, qArr] = iterate_nhmm(u, k, a, b, x0, q0, N);

    % 只取后 keep 步
    xPlot = xArr(end-keep+1 : end);
    qPlot = qArr(end-keep+1 : end);

    % 在 (2×2) 的网格中第 i 个子图
    subplot(1,3,i);
    plot(xPlot, qPlot, '.', 'MarkerSize', 3);
    xlabel('x'); ylabel('q');
    title(['b = ', num2str(b)]);
    grid on; box on;

    % xlim([-20,20]); ylim([-1,0.5]); end
end

sgtitle('不同 b 值下的 (x,q) 吸引子');
